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CO ! ABSTRACT 
O ■ 

' Anomalous microwave emission (AME) has been observed by numerous experiments in the frequency range ~ 10-60 GHz. Using Planck maps 
* and multi-frequency ancillary data, we have constructed spectra for two known AME regions: the Perseus and p Ophiuchi molecular clouds. The 
>^ spectra are well fitted by a combination of free-free radiation, cosmic microwave background, thermal dust, and electric dipole radiation from 
small spinning dust grains. The spinning dust spectra are the most precisely measured to date, and show the high frequency side clearly for the 
I ' first time. The spectra have a peak in the range 20-40 GHz and are detected at high significances of IV.lcr for Perseus and 8.4cr for p Ophiuchi. 
. . I In Perseus, spinning dust in the dense molecular gas can account for most of the AME; the low density atomic gas appears to play a minor role. 
^ . In p Ophiuchi, the ~ 30 GHz peak is dominated by dense molecular gas, but there is an indication of an extended tail at frequencies 50-100 GHz, 
which can be accounted for by irradiated low density atomic gas. The dust parameters are consistent with those derived from other measurements. 
We have also searched the Planck map at 28.5 GHz for candidate AME regions, by subtracting a simple model of the synchrotron, free-free, and 
thermal dust. We present spectra for two of the candidates; S140 and S235 are bright H ii regions that show evidence for AME, and are well fitted 
by spinning dust models. 

Key words. ISM: general - Galaxy: general - Radiation mechanisms: general - Radio continuum: ISM - Submillimeter: ISM 
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1. Introduction 

Anomalous microwave emission (AME) is an additional 
component of diffuse foreground emission that cannot be 
easily explained by synchrotron, free-free, or thermal dust 
emission. AME has been observed by numerous experiments 
over the frequency range -10-60 GHz and is found to be very 
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closely correlated with far infrared (FIR) emis sion associ- 
ated wi th thermal e niissio n from dust grains (Kogut et al.J 
199 d iLeitgh et alJ |1997|; Ide Qliveira-Costa et al." 1997^: 
jLag 



Bandav et al 



lay .-^ - ^ - ^ 

2004: 'Finkbeiner"2QQ4:"Finkbeineretal.' '2004': 'Davies et al 
2006; Dobler & Finkbeiner 2008; Miville-Deschenes etal 
2008; Gold et al.1 120111: lYsard et al.1 |2010). Electric dipole 
radiation from small rapidly spinning dust grains, or "spinning 
dust," is thought to be emitted in the microwave region of 
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the spectrurrQ. Theoretical models predict a peaked spec- 
trum, at a range of possible frequencies (-10-150 GHz) 
depending on the properties of the dust gr ains and th eir envi - 
ronment fc raine &Lazarian 1998; AH-H annoud et al.l l2009l: 
iDobler et all .2009: .Ysard & Verstraete 201Qb . The spectrum 
is very different from those of the traditional diffuse emission 
expected at these frequency ranges; for example, free-free, 
synchrotron, and thermal dust emission have power-law-like 
spectra at microwave frequencies. For these reasons, spinning 
dust emission has become the best explanation for the AME, 
although other physical mechanisms could still be contributing 
at s ome level, such as hot free-free em ission (iLeitch et al.1 
1199 7), hard svnchrotron radiation (Bennett et al. I 12003 ). or 
magneto-dipole emission dPraine & Lazarianll 19991) . 

The most direct evidence for spinning dust grains has come 
from dedicated observations of known dust clouds in the in- 



terstellar medium (Einkbeiner et al. 2002; Casas sus et al.l 12006 : 
Dickinson etaD[2009: Scaife et al. 2009: Dickinson et al.ll2009 . 



20101) . Per haps the best exampl e is that of the Perseus molec 
ular cloud. Watson et al. ' (^2005^ observed the region with the 
COSMOSOMAS telescope covering 1 1-17 GHz. By combining 
data from radio surveys, WMAP, and DIRBE, they were able to 
show very clearly a strong excess at frequencies -10-60 GHz 
and a shape similar to spinning dust models. Another clear ex- 
ample is that of the photodissociation region (PDR) associated 
with the p Ophiuchi molecular dust cloud. Casassus et al. ( 2008) 
observed the region centred at (/, b) = (353?0, -hl7?0) at 31 GHz 
with the Cosmic Background Imager (CBI) and found a close 
connection between the centimetre-wavelength emission and 
the dust emission. Eurther analysis of the spectrum confirmed 
that there was significant excess emission that was consistent 
with spinning dust models. AME has also been d etected within 
a star -forming region in the galaxy NGC6946 ('Murp hy et al.l 
120101). which ap pears to be well-fitted by spinning dust models 
dScaife et al. 2010). 

PlancJ^ data now represent the next step in measuring and 
understanding spinning dust emission. The additional data at 
28.5, 44.1, and 70.3 GHz (where spinning dust is expected to 
emit strongly) and the high frequency data at 100-857 GHz al- 
low an accurate model for the thermal emission to be removed. 
This is crucial in defining a precise spectrum for the spinning 
dust, particularly at the higher frequencies where the smallest 
dust grains dominate. Accurate spectra for spinning dust could 
potentially prove to be an important diagnostic tool for measur- 
ing the properties of the dust (e.g., the size distribution, density 
and average electric dipole moment) as well as the physical envi- 
ronment in the surrounding region (e.g., the interstellar radiation 
field). 

In this paper, we study two known AME regions, the Perseus 
and p Ophiuchi molecular cloud regions, and provide spectra 
over the frequency range 0.4-3000 GHz. We show, for the first 
time, an accurate residual spectrum after removal of the free- 
free, cosmic microwave background (CMB), and thermal dust 
components. The AME is clearly detected as a convex spectrum 
that is well-fitted by theoretical models of spinning dust emis- 
sion. We investigate the parameters for the spinning dust to see 



^ lEricksonI (11957 ^ was the first to suggest the possibility of radio 
emission from spinning dust grains. 

^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 



whether they are reasonable, given the conditions in these envi- 
ronments. The Planck data, in combination with ancillary data, 
are also used to search for new regions of AME. We present two 
such regions and show that these newly discovered regions can 
also be fitted with a contribution from spinning dust grains. 

In this paper, Sects.O and [3] describe the Planck and ancil- 
lary data that were used to measure the spectra of AME regions. 
Section |4l presents the results for the Perseus molecular cloud 
and Sect.[5]for the p Ophiuchi cloud. Section[6]describes prelim- 
inary modelling of the AME in terms of spinning dust. Section[7] 
presents the search for new AME regions and gives two clear 
examples. Conclusions are given in Sect.[8l 



2. Planck data 



Planck (' Tauber et al.ll2010l: iPlanck Collaboration! l20 1 la*) is the 
third-generation space mission to measure the anisotropy of the 
cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular r esolution from 31' to 5'. The Low Erequency 
Instrument (LEI: iMandolesi et al.ll20TQl: iBersanelh et aLlboiol: 
iMennella et al.ll201 ih covers the 30. 44, and 70 GHz bands with 
amplifier s cooled to 20 K. The High Erequenc y Instrument (HEI; 
Lamarr e et al ] 120101: iPlanck HEI Core TeamI [201 la) covers the 
100, 143, 217, 353, 545, and 857 GHz bands with bolome- 
ters cooled to O.I K. Polarisation is measured in all but the 
highest two bands (iLeahv et al.l 120101: iRosset et al] l2010l) . A 
combination of radiative cooling and three mechanical cool- 
ers produces the temperatures ne eded for the detectors and op- 
tics (iPlanck Collaborationll201 Ibl) . Two data processing centres 
(DPCs) check and calibrate the data and make maps of the sky 
(IPlanck HEI Core Team 2011b; Zaccheietal. 2011). Planck's 
sensitivity, angular resolution, and frequency coverage make it a 
powerful instrument for Galactic and extragalactic astrophysics 
as well as cosmology. Early astrophysics results are given in 
Planck Collaboration, 201 Ih-z. 

In this paper, we start with the Planck (non CMB -subtracted) 
maps at nominal frequencies of 28.5, 40.1, 70.3, 100, 143, 217, 
353, 545, and 857 GHz (see Table [T]). For the spectral analyses, 
we smoothed the HEALPix maps to a common angular reso - 
lution of r and degraded to A^side = 512 ( Gorski et al. 2005"). 
We assu me Gaussian beams and use th e average beamwidth s 
given in IPlanck HEI Core TeamI (l2011bh : IZacchei et al.l (l201ll) . 
Details of the beam, such as the exact shape and variation of 
the beam across the sky, do not significantly aflTect the results 
presented here, due to the applied smoothing which dominates 
the spatial response of the maps. Furthermore, extended emis- 
sion is less sensitive to the details of the beam. We convert 
from CMB thermodynamic units to Rayleigh- Jeans brightness 
ter nperature units u sing the standard conversion factors given 
in IPlanck HEI Core Team! (1201 IbD and IZacchei et all (1201 ID . 
Colour corrections to account for the finite bandpass at each fre- 
quency are applied during the modelling of the spectra for each 
object; these are typically a few per cent for LEI and 10-15 % 
for HEI. The 100 GHz data are significantly contaminated by the 
CO / = 1 ^ line at 1 15 GHz. We apply a no minal correcti on to 
the 100 GHz map using the CO survey of Dame et al.l (1200 ll) . We 
multiply the integrated 1° -smoothed line intensity map (in units 
of K km s-i) by 14.2 ('Planck HEI Core Teamll2QTlbl) and sub- 
tract it from the 100 GHz map. Instead of attempting to correct 
for other lines that may contaminate the 217 and 353 GHz data, 
we do not include these channels in the modelling of spectra if 
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Table 1. Summary of the data. 



Frequency 


Telescope/ 


Angular resolution 


Coverage 


Notes 


[GHz] 


survey 


[arcmin] 






0.408 


JB/Eff/Parkes^ 


^ 60 


Full- sky 


NCSA desourced/destriped version ^ 


0.820 


Dwingeloo'' 


72 


Northern sky (-7° <6 < +85°) 


Available on the web^ 


1.420 


StockertA^illa-Elisa^ 


36 


Full- sky 


Courtesy of W. Reich 


2.326 


HartRAO^ 


20 


Southern sky (-83° <6 < +13/32°) 


Courtesy of J. Jonas (priv. comm.) 


11-17 


COSMOSOMAS^ 


^ 60 


Northern sky (+24° <6 < +44°) 


Filtered on large angular scales 


22.8 


WMAP 7-yr^ 


^ 49 


Full- sky 


1 ° smoothed version^ 


28.5 


Planck LFr 


32.'65 


Full- sky 




33.0 


WMAP 7-yr^ 


^ 40 


Full- sky 


1 ° smoothed version^ 


40.7 


WMAP 7-yr^ 


^ 31 


Full- sky 


1° smoothed version^ 


44.1 


Planck LFr 


27.^92 


Full- sky 




60.7 


WMAP 7-yr^ 


^ 21 


Full- sky 


1° smoothed version^ 


70.3 


Planck LFr 


13.'01 


Full- sky 




93.5 


WMAP 7-yr^ 


^ 13 


Full- sky 


1° smoothed version^ 


100 


Planck HFF 


9.37 


Full- sky 




143 


rtl 7 X TT7T / 

Planck HFF 


7. '04 


Full- sky 




217 


Planck HFF 


4.'68 


Full- sky 




353 


Planck HFF 


4.'43 


Full- sky 




545 


Planck HFF 


3.'80 


Full- sky 




857 


Planck HFF 


3.'67 


Full- sky 




1249 


COBE-DIRBE'' 


;^40 


Full- sky 


LAMBDA website^ 


2141 


COBE-DIRBE'' 


=^40 


Full- sky 


LAMBDA website^ 


2997 


COBE-DIRBE'' 


^40 


Full- sky 


LAMBDA website^ 



Notes . ^^^lHaslametal.1 (119821) http : //lambda^sf c . nasa.gov/ IB erkhuiiseni (119721) ht tp : //ww w . mpi f r - bonn . mp g . de/ survey . html , ^^^ lRej 

I Jonasetal.1 (119981) IWatson et al.l (l2005l) [Jarosik et all (I20T1I) IZacchei et al.l 



1982 ): Reich & Reich (1986); Reich etal. (200l| 



2011) ^^\Planck HFI Core Team (2011b.) ^^\Hauser et al 



(119981) 



they are observed to be in excess by more than 2cr, relative to the 
best-fit model when omitting these frequencies. Uncertainties 
are taken to be the nominal uncertainties in absolute calibration 
and gain (given in the aforementioned papers), which are added 
in quadrature. If these are below 3 %, we conservatively assume 
an overall uncertainty of 3 %. 

Jack-knife tests, in which the data were split into two 
halves (in terms of time and also HFI detectors; see 
iPlanck HFI Core Teamll201 Ibl) . provided a check of the consis- 
tency of the Planck data. The tests confirmed that the results 
were consistent with each other to within a small fraction of the 
derived uncertainties. 



(Paradis et al., in prep.). We note that there are significant base- 
line uncertainties (i.e., off'sets) in these maps, but these will not 
aff'ect our results since we are subtracting a local background for 
each object, which removes any such off^set. Th e data include 
the full-sky map at 4 08 MHz dHaslam et al.lll982h . the 820 MHz 
Dwingeloo sur vev dBe rkhuiisen 1972^, the new full- sky map 
at 1420 MHz ([ReichP l982r Reich & Reich 1986; Reich et al' 
2001), and the 2326 MHz HartRAO southern survey (iJonas et af 



.1998) . The 1420 MHz map was multipHed by a factor of 1.55 to 
account for the main beam to full beam ratio (iReich & ReichI 
1988), and to bring the calibration in line with the other data. 
We assume a 10 % overall calibration uncertainty for these data. 



3. Ancillary data 

Although Planck observes over a wide frequency range (30- 
857 GHz), it is still important to include ancillary data to extend 
the frequency range into the FIR, and in particular, to lower fre- 
quencies. To supplement Planck data we therefore use a range 
of publicly available surveys and data available within the team. 
Table[T] summarises the data used in the analysis and the origin of 
the maps that were used. In the following sections, we describe 
the various data sets in more detail. 

3.1. Low frequency radio data 

The radio data consist of a number of low frequency surveys in 
the range 0.408-2.326 GHz (Table [B. We used the most up-to- 
date versions of the data sets, and where necessary, regridded 
the maps into the HEALPix format (iGorski et al.ll2QQ5[) by using 
a procedure which computes the surface intersection between 
individual pixels of the survey with the intersecting HEALPix 
pixels. This procedure has been shown to conserve photometry 



3.2. COSMOSOMAS data 

The COSMOSOMAS experiments ("Gallegos etaP 
l200lh consist ed of two circu lar s canning instruments 
COSMO 11 (lHildebr andtetal.1 l2QQ7h and COSMO 15 
(iFernandez-Cerezo et al.i i2006) which operated at the Teide 
observatory (altitude 2400 m, Tenerife). They have produced 
0?8-l?l resolution maps of 10,000 deg^ in four frequency 
bands (10.9, 12.7, 14.7, 16.3 GHz) that are ideal for filling 
the gap between low frequency radio surveys and the lowest 
frequency channels of WMAP and Planck. The average sensi- 
tivity of each map is in the range 80-120 m^K per beam,^nd 
the observing strategy is described in iGallegos et al.l (l2Q0ll) . 
As COSMOSOMAS was a ground-based experiment, there is 
significant contamination from atmospheric and receiver l/f 
noise, which needs to be filtered out. This results in a loss of 
response to angular scales larger than 5° and introduces low 
level correlations over the scan ring. This must be taken into 
account when considering diff'use emission at angular scales 



3 



Planck Collaboration: Planck Early Results. XX. New light on anomalous emission from spinning dust grains 



> 1° and this will be described in SectS A 10% overall 
calibration uncertainty is assumed. 

3.3. WMAPdata 



WMA P 7-year data are included in the analysis (iJarosik et alJ 
l2QTTh . The data span 23-94 GHz and thus complement Planck 
data, particularly the K-band (22.8 GHz) channel. We use the 1°- 
smoothed maps available from the LAMBDA website. We apply 
colour correctio ns to the central frequencies using the recipe de- 
scribed by Ja rosik et al.l (|2003) and assume a conservative 3 % 
overall calibration uncertainty. 

3.4. DIRBE data 

To sample the peak of the blackbody curve for tempera- 
tures > 15K we include the COBE-DIRBE data at 240 ^tm 
(1249 GHz), 140 ^tm (2141 GHz) and 100 ^tm (2997 GHz). The 
DIRB E data are the Zodi- Subtracted Mission Average (ZSMA) 
maps (iHauser et al.lll998l) regridded into the HEALPix format. 
Colour corrections were applied as described in the DIRBE ex- 
planatory supplement version 2.3. Data at higher frequencies 
are not considered, since these will be dominated by transiently 
heated grains which are not in thermal equilibrium with the in- 
terstellar radiation field and therefore cannot be modelled with a 
single modified blackbody curve. Furthermore, at wavelengths < 
40 |Ltm the spectrum is dominated by many emission/absorption 
lines. 



4. Perseus Molecular Cloud 

4.1. Introduction and Maps 

The Perseus molecular cloud complex is a relatively nearby gi- 
ant molecular cloud of 1.3 x 10"^ M© at a distance of 260 pc 
(ICernicharo et al.|[T985l) having an angular extent of 6° x 2°. It 
contains a chain of six dense cores (B 5, IC 348, B 1, NGC 1333, 
L 1455 and L 1448) with Ay > 2 mag running northeast- 
southwest and with two main centres of medium-mass star 
formation associated with the reflection nebulae IC348 and 
NGC 1333 near each end (see Fig. [T]). As a nearby example of 
low-to-medium mass star formation it has been extensively stud- 
ied, with the Spitzer c2d Lega cy Progra m (Evans et al. 2003) 
and the COMPLETE survey ( Ridge et al.ll2Q03l) . which provide 
a wealth of infrared and spectral line data to understand condi- 
tions throughout the cloud. 

Anomalous emission was first detected as an excess in this 
regio n of the sky via the Tenerife experiment (iDavies et al.l 
1 19871) at 5° angular scales, a nd it contributed most of th e 
dust-correlated signal found by Ide Qliveira-Costa et al.l (Il999l) . 
With the 1° beam and 11-1 7 GHz frequency b ands of the 
COSMOSOMAS experiment (Gallegos et al.l l200lh it was pos- 
sible to locate and identify the Perseus molecular cloud as the 
strong source of this emission (IWatson et al.ll2005|). Follow- up 
observations with the VSA interferometer fTibb s et al.ll20TQh at 
10' resolution detected five "hot-spots" of anomalous emission, 
one of which corresponds to the dense region (Ay > 10) around 
the IC 348 core. The other four hot-spots are associated with 
nearby dense dust knots that lie on a warm ring of dust emission 
(G159.6-18.5), about 90' to the south-east of IC 348. This ring, 
which is seen at 100 |Lim, is thought t o be warmed by the OB 
star HD 278942 (ide Zeeuw et al .11 19991) . which is also responsi- 
ble for a weak Hii region in the centre of the ring, and which 



can be seen in the low-frequency continuum surveys with an in- 
tegrated flux density of a few janskys. The AME appears to have 
littl e or no polarisation with an upper limit of 1.8 % at 33 GHz 
dLopez-Caraballo et aLllMoh . 

Planck maps of the Perseus molecular cloud region, cover- 
ing 30-857 GHz, are shown in Fig. (TJ along with the 1.4 GHz 
and Ha maps. For display purpose s the Planck maps have 
been CMB- subtracted a s descri bed in IZacchei etal.l (|201ll) and 
iPlanck HFI Core TeamI (1201 Ibh and the units are thermody- 
namic kelvins. The strong dust-correlated AME at 30-70 GHz 
is evident; it has no obvious counterpart at 1 .4 GHz but corre- 
lates well with the higher (> 100 GHz) Planck frequencies, which 
are dominated by thermal dust. The majority of the AME comes 
from the northeast end of the molecular cloud, which harbours 
the dense IC 348 reflection nebula. Some weak free-free emis- 
sion is seen to the east and the contribution from the extragalactic 
radio source 4C-h32.14 is clear and must be removed (Sect. l4.2l) . 
The Hof (and 1 .4 GHz emission) morphology is not co-located 
with the emission seen at 30-70 GHz, but is likely to be signifi- 
cantly absorbed by the intervening dust. Nevertheless, t he bright- 
est Hof pixels predict ^ 0.25 mK at 30 GHz (Dickin son et al.l 
12003) . a factor of ten below the AME at 30 GHz. It is interest- 
ing that there is considerable bright structure observed in the 
Planck HFI maps, which is not always associated with strong 
AME. This is likely to be a consequence of a diff'erent environ- 
ment and dust grain population. 



4.2. Spectrum 

To reduce correlated 1 // noise from the atmosphere/receiver the 
COSMOSOMAS time-ordered data have been filtered with the 
suppression of the first seven harmonics in the EFT of the circu- 
lar scans. To compare with flux densities from other instruments 
one must either correct for the flux lost in filtering or filter all 
the other data in the same way. With extended sources such as 
Perseus there is uncertainty in estimating a single correction fac- 
tor to take into account the filtering process, which in turn di- 
rectly aff'ects the shape of the frequency spectrum. Filtering all 
the data in the same way retains the spectral shape and ensures 
that the spatial frequencies used to form the spectrum are com- 
mon between the diff'erent data sets. 

The first step is to smooth all the data to the same resolu- 
tion as the lowest frequency channel of the COSMOSOMAS 
data, i.e., a FWHM of 1?12. Each data set is then fil- 
tered using the COSMOSOMAS reduction pipeline and binned 
into a map as was done with the original COSMOSOMA S 
data (iFernandez-Cerezo et al.l l2006l: iHildebrandt et al.1 l2007h . 
Example maps at 10.9, 28.5, 70.3, and 857 GHz are shown in 
Fig. [21 To estimate the flux density at each frequency the filtered 
maps were compared to a filtered model map, which was passed 
through the same scan strategy. For the model we use an ellip- 
tical Gaussian (1?6 x 1?0 with major-axis PA 51°), centred at 
(RA,D ec) = (55?33 +31?6 7) or (/,/?) = (160?26, -18?62), fol- 
lowing l^tsonetal] (l2005l) : this gives a solid angle D = 5.52 x 
10-4 sr. We win refer to this as AME-G160.26-18.62. There is 
a flat-spectrum radio quasar, 4C-h32.14 (NRAO 140), within 1°, 
with a flux density of ^ 1-3 Jy, which is seen in the Planck Early 
Release Compact Source Catalogue (ERC SC) at frequencies up 
to 143 GHz (P lanck CoUaborationll20TT^ . We use the flux den- 
sities and beam widths taken from the ERCSC and subtract the 
two-dimensional Gaussian profile from the maps before smooth- 
ing and filtering. Meanwhile in the COSMOSOMAS maps we 
must extrapolate the flux density of 4C-h32.14 to the appropriate 
frequencies and simulate its contribution to the filtered map; we 
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Fig. 1. Maps of the Perseus molecular cloud region at their original angular resolution. From left to right, top row: Planck 28.5; 
44.1; 70.3 and 100 GHz. Bottom row: Planck 143 and 857 GHz; 1.4 GHz; and Ha. The maps cover 5° x 5° centred on (/,/?) = 
(160.°26, -18.°62) and have linear colour scales. The graticule has 1° spacing in Galactic coordinates. The FWHM of the elliptical 
Gaussian model used to fit the flux density in the filtered maps (see text) is shown. The strong AME is evident at 30-70 GHz. 




Fig. 2. 16° X 16° maps after filtering with the COSMOSOMAS data reduction pipeline. The maps are centred at (/, b) = (160°, - 15°) 
and the graticule spacing is 2° in Galactic coordinates. Erom left to right: COSMOSOMAS 10.9 GHz; Planck 28.5 GHz; 70.3 GHz; 
and 857 GHz. The elliptical Gaussian model used to calculate the spectrum is indicated. 



assume 2.0 Jy at 11-19 GH z, which is in good agreement with 
NVSS (ICondon et al.lll998l) and recent 2.7 GHz measurements 
made with the Eff'elsberg telescope (iReich & Reich"'2009). The 
contribution of 4C+32.14 is 5 % of the integrated flux density, 
so the final results are not very sensitive to the exact treatment of 
the subtraction or the variability of the source. Both the smooth- 
ing and filtering result in correlated noi se betwee n the pixels in 
the map. A correlation analysis (e.g., D avies et al .11200 6) is used 
to calculate the flux density taking into account the correlated 
noise. The co variance matrix is estimated from the size of the 



smoothing kernel and the point source response to the scan strat- 
egy. A simulated map is produced for the source, with an inte- 
grated flux density of 1 Jy, and a maximum likelihood method 
is used to find the amplitude and error for the given frequency 
channel data within a 19 x 19 pixel box (6?3 x 6?3) around the 
source position. 

As a control, we tested the simulation and filtering pipeline 
by applying the technique to the extended California nebula H ii 
region (NGC 1499), approximately 6° to the north of Perseus 
at (/,/?) = (160?60,-12?05), as seen in Eig. [2l The emission 
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Fig. 3. Spectrum of the California nebula (NGC1499), measured 
using the filtered flux method (see text). A simple power-law fit 
is shown. 

at radio/microwave wavelengths is dominated by optically thin 
free-free emission. Flux densities at frequencies >100 GHz were 
afl'ected by dust emission lying adjacent to the north-west of 
NGC 1499 and thus were not included. The spectrum is shown 
in Fig. O A simple power-law fit, with a best-fitting spec- 
tral index a = -0.05 ± 0.01, is overplotted. The fit is rel- 
atively good (x^/dof= 1.23) but there is some evidence for 
AME at frequencies ^ 30-100 GHz. Furthermore, the spectral 
index is slightly flatter than the value expected for optically 
thin free- free emission over this frequency range, a ^ -0.1 
to -0.15 (iDickinson et al.l l2003h . Indeed, AM E has been ob- 
served to be significant in some Hii regions ( Dickinson et al.l 
120091: Pibdorovic et"ani2010h . We confirmed the spectral shape 
with a standard aperture photometry analysis (as described in 
Sect.O and will investigate this in more detail in a future publi- 
cation. Note that the excess at 30-100 GHz disappears when the 
spectrum is calculated away from the nearby dust feature to the 
north-west of NGC 1499. Given the consistency of the flux den- 
sities between the frequencies, and that we have obtained similar 
results using standard aperture photometry, we are confident that 
the spectrum calculated using the filtered maps is robust. 

The spectrum for AME-G160.26- 18.62 in Perseus is shown 
in Fig. 131 The flux densities and associated errors are listed 
in Table O The spectrum is well sampled across the radio, 
microwave and FIR regi mes and i s a significant improvement 
on that presented in Watson et alj (12005), with the additional 
Planck data allowing a much more accurate spectrum of the 
AME to be extracted. The low-frequency (< 2 GHz) data show 
a flat spectrum consistent with free-free emission while the high 
frequencies (> 100 GHz) are dominated by thermal dust emis- 
sion. The excess at ^ 10-70 GHz is evident and has a peaked 
(convex) spectrum, with a maximum at 25 GHz. 

The data allow us to fit a multi-component parametric model 
to the flux density spectrum. The model consists of four compo- 
nents: free-free emission; thermal dust emission; a CMB fluctu- 
ation; and spinning dust emission. The sum is 



■ ^ ff + ^ td + ^ cmb + ^ 



' sp- 



(1) 

The free-free flux density, S^, is calculated from the brightness 
temperature, T^, based on the optical depth, tq, using the stan- 
dard formulae: 

_2kT^ 
^ff o ' 
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Fig. 4. Spectrum of AME-G160.26-18.62in the Perseus molec- 
ular cloud. The best-fitting model consisting of free-free (orange 
dashed line), spinning dust, and thermal dust (light blue dashed 
line) is shown. The two-component spinning dust model consists 
of high density molecular gas (magenta dot-dashed line) and low 
density atomic gas (green dotted line). 
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Fig. 5. Spectrum of AME-G160.26-18.62in the Perseus molec- 
ular cloud after subtracting the best-fit free-free, CMB and ther- 
mal dust components. A theoretical spinning dust model, con- 
sisting of two components is shown as the black solid line; the 
magenta dot-dashed line is for high density molecular gas and 
the green dotted line is for low density atomic gas, which is a 
very small contribution in this case (see Sect.©. 



where k is the Boltzmann constant, Q is the solid angle, and v is 
the frequency, 

Tff = Te(l - e-'% (3) 
and the optical depth, Tff , is given by 

Tff = 3.014 X 10"^r;^-^y-^EMgff, 



(4) 



where Tg is the electron temperature (in units of K), EM is the 
emission measure (in units of cm"^ pc) and is the Gaunt fac- 
tor, which is approximated as 



^ /4.955xl0-2\ ^ 

gs = M — +1-5 HTe). 

y/GHz 



(5) 
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Table 2. Flux densities for AME-G160.26- 18.62 in the Perseus 
molecular cloud and residual fluxes when free-free, CMB and 
thermal dust components are removed. The fitted spinning dust 
model consists of two components. 



Frequency 


Flux density 


Flux density residual 


[GHz] 


[Jy] 


[Jy] 


0.408 


9.7 + 3.3 


0.2 + 3.4 


0.82 


9.4 + 4.7 


0.5 + 4.8 


1.42 


8.0+ 1.7 


-0.5 + 1.8 


10.9 


16.1 + 1.8 


9.4 + 2.0 


12.7 


20.0 + 2.2 


13.3 + 2.4 


14.7 


28.4 + 3.1 


21.8 + 3.3 


16.3 


35.8 + 4.0 


29.3 + 4.2 


22.8 


39.8 + 4.2 


33.4 + 4.5 


28.5 


38.1 +4.0 


31.5 + 4.4 


33.0 


36.7 + 3.9 


29.7 + 4.4 


40.9 


30.8 + 3.3 


22.6 + 4.0 


44.1 


28.5 + 3.2 


19.6 + 4.0 


61.3 


27.4 + 3.4 


10.9 + 6.4 


70.3 


32.2 + 3.9 


8.9 + 9.0 


93.8 


63.0 + 7.8 


7 + 21 


100 


78 + 15 


10 + 29 


143 


202 + 22 


-25 + 80 


217 


1050 + 130 


120 + 320 


353 


3860 + 470 


-600 + 1400 


545 


15300 + 2100 


-800 + 4900 


857 


48700 + 6100 


-1000+ 14000 


1250 


93000 + 13000 


-2400 + 28000 


2143 


117000 + 15000 


7000 + 35000 


2997 


53600 + 6700 


-2000 + 20000 



100 GHz and is not included in the fit. Similarly, the 217 GHz 
value is not included in the fit, due to the possible small con- 
tamination from the CO / = 2 ^ 1 line. From the initial spec- 
tral fit the spectral index is found at each frequency, which is 
then used to calculate the appropriate colour corrections for each 
band in an iterative manner. After applying the colour correc- 
tions (typically a few per cent correction), the spectrum is fit- 
ted again to obtain the final flux densities, which are given in 
Table O The uncertainty in the fitted model for the spectrum is 
computed by propagating the errors of each of the parameters 
and combining them with the other sources of error, i.e., back- 
ground/noise residuals and an overall calibration error. The fit is 
very good with;^^/dof = 0.95. The best-fitting parameters are: 
ElVl = 148 ± 13.5 cm-6pc for free-free; = 18.5 ± 0.6 K and 
= 1.65 ± 0.08 for the thermal dust; and a neghgible CIMB 
contribution of ATqmb = -6 + 66 |LtK. 

In Fig. [5] we show the residual spectrum after subtracting the 
free-free, CIMB and thermal dust models. The additional uncer- 
tainty from the subtraction of the models is added in quadrature 
to the flux density errors to find the final uncertainties on the 
residual flux densities, which are listed in Table [2j This repre- 
sents the most precise spectrum of AIME to date. The weighted 
average over the range 10 to 94 GHz gives a detection signifi- 
cance level of 17.1 (T. The peak of the spinning dust component 
is centred at ^ 25 GHz and the low and high frequency sides are 
well-defined in the range 10-90 GHz. A physical model for the 
spinning dust is presented in Sect.O It consists of two compo- 
nents (atomic and molecular), which are overplotted in Fig. \5\ 
and is seen to be an excellent fit to the data. 



It is the Gaunt factor that results in a slight steepening of the 
free-free spectrum with frequency, which is particularly pro- 
nounced at > 100 GHz. The thermal dust is modelled as a single- 
component, modified blackbody curve, y^B(y, Td), which we 
normalise using an optical depth at 250 |Ltm (1.2 THz), T250, as 

= 2''5^^wI^^25o(v/1.2THz)/'Q, (6) 

where h is the Planck constant and /S is the emissivity index. The 
CMB anisotropy is evaluated as the diff'erential of the black- 
body temperature function, in thermodynamic units (ATqub), 
converted to flux density units, 

^cmb - I — -2 — I ^^CMB- (7) 

The spinning dust model is based on an assumed theoretical 
emissivity curve, jy (in units of Jy sr"^ cm^ H"^), normalised by 
the average column density per H nucleon, A^h, and solid angle, 

a, 

S,^=Nnjyn. (8) 

As described in Sect .[6l we use the SPDUST code 
i-Haimoud et al ]|200i) to model the spinning dust spectrum. 
The best-fitting model is plotted in Fig.|4l The fitted param- 
eters are EM for free-free, diff'erential temperature of CMB per- 
turbation (ArcMB), dust temperature (T^), dust emissivity index 
(j3) and 250 jam opacity (T250). The spinning dust model con- 
sists of two components, to represent dust associated with dense 
molecular gas and low density atomic gas, and will be discussed 
in Sect.[6j The electron temperature is assumed to have a typi- 
cal value of 8000 K for gas in the solar neighbourhood. A con- 
servative 3 % uncertainty in the CO correction is assumed at 



5. p Ophiuchi Molecular Cloud 

5.1. Introduction and Maps 

The p Ophiuchi molecular cloud (e.g., lEncrena 
iKulesa et al.ll2005l: lYoung et al.ll20Q6l) Hes in the Gould Belt of 
the closest molecular complexes, at a distance D = 135+15 pc 
(Ferryma n et al.|[T997h . It is undergoing intermediate-mass star 
formation. Ultra-violet radiation from its hottest young stars 
heats and dissociates exposed layers, but does not ionise bulk hy- 
drogen. The most prominent FDR in p Ophiuchi is the p Oph W 
filament, with nu = 10"^ - 10^ cm~^, which h as been studied by 
iLiseau et al.l (|1999) and Habart et alJ (l2003h . p Oph W is ex- 
cited by HD 147889 (a B2,3 IV binary), the earhest star in the 
p Oph star formation complex. The line of sight to HD 147889 
can be taken to characterise the physical conditions in the ambi- 
ent diff'use cloud, with A^h ~ 2 x 10^^ c m~^, fin ~ 400 cm~^; 
no ac curate values are available for A^H2 ( van Dis hoeck & Blackl 
119891: [Kazmie rczak et al. 2010). The bulk of the mass in the 
p Ophiuchi cloud is situate d in the Oph A molecular core (e.g., 
IWootten et all I1978L IT979b . located roughly 10' southwest of 
p Oph W. 

Only faint radiation from the Ray leigh- Jeans tail of ~ 10- 
100 K dust is expected at wavelengths long wards of ~3 mm. Yet 
CBI observations revealed that the p O ph W PDR is surpris - 
ingly bright at centimetre wavelengths ( Casassus et al.l [2008h . 
The WMAP 33 GHz image confirms that the general location 
of the signal seen by CBI is off'set by 10' from the centroid of 
the WMAP 94 GHz data, located on the Oph A core. The sig- 
nal seen at centimetre wavelengths has no 5 GHz counterpart; 
instead the optically thin free-free emission is located in the im- 
mediate vicinity of HD 147889. The other early-type stars in the 
complex are S 1 (B4 V) and SR 3 (B6 V), and each is surrounded 
by bright IR emission peaks, which lack any radio counterpart. 
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Frequency [GHz] 

Fig. 7. Spectrum of the Orion nebula (M42) Hii region using 
aperture photometry. The model consists of optically thin free- 
free emission (orange dashed line) and thermal dust emission 
(light blue dashed line). The inner panel shows the residual spec- 
trum after removal of the model, indicating the consistency be- 
tween Planck and WMAP data. 



Planck maps of the p Ophiuchi molecular cloud region cov- 
ering 30-857 GHz, along with the 1.4 GHz and Ha maps, are 
shown in Fig. [6l Strong AME is evident at 28.5-44.1 GHz, 
with the bulk of the emission falling within 0.°5 of the cen- 
tral position. Eor display purposes the Planck maps have 
been CMB- subtracted a s descri bed in lZacchei et al.l (1201 1) and 
iPlanck HFI Core TeamI (l2011bh and the units are thermody- 
namic kelvins. The H ii region to the west is also clearly detected 
in the Planck LEI maps, but is well separated from the bulk of 
the AME. In the location of the strongest AME, there is a val- 
ley in intensity seen at low frequencies, supported by little Ha 
emission. The AME is well correlated with Planck HEI maps, 
indicating an association with dust. 

5.2. Spectrum 

The p Ophiuchi molecular cloud (at declination -24°) lies out- 
side the range of the COSMOSOMAS survey. We therefore 
chose to use simple aperture photometry on the unfiltered maps 
to calculate the spectrum. The maps were first smoothed to a 
common resolution of 1° and converted to units of Jy pixel" ^ 
The brightness in each pixel was then summed over a given aper- 
ture to give the flux density. We chose a circular aperture centred 
at (/, b) = (3 5 3? 05, 16?90) with a radius of 60'; this corresponds 
to a solid angle Q = 9.57 x 10^ sr. We will refer to this source 
as AME-G353.05-hl6.90. An estimate of the background bright- 
ness is subtracted using the median value in a circular annulus 
between an inner radius of 80' and an outer radius of 120', as 
indicated in Eig.O This removes local background emission and 
any residual ofl'sets in the maps. Errors are estimated from the 
standard deviation in the background annulus, which was ver- 
ified with simulations of random sources of known flux den- 
sity that were injected into the sky maps. The results were not 
strongly sensitive to the exact choice of background annulus; 
all values were within the Icr uncertainties. This indicates that 
the emission in the aperture is bright relative to the local back- 
ground fluctuations. We also compared two halves of the data as 



a jack-knife test. Results were found to be consistent to within a 
fraction of the derived uncertainties. 

We tested the aperture photometry by measuring the spec- 
trum of well-known, bright objects. We chose bright H ii regions, 
which are expected to be dominated by thermal bremsstrahlung 
(free-free) radiation at radio frequencies and thermal dust in the 
sub-mm/EIR. Eig. [7] shows the spectrum of the Orion nebula 
(M42) with an aperture of radius 60' and background annulus 
between radii of 80' and 120'. At frequencies below a few gi- 
gahertz, the free-free emission is becoming optically thick, with 
a turnover frequency of ^ 2 GHz. At frequencies -10-60 GHz, 
the spectrum is well-fitted by a simple optically thin free-free 
spectrum (a ^ -0.1). The WMAP and Planck data points are 
in good agreement to within ^ 1 %; the residuals are consistent 
with zero and within a fraction of the assumed uncertainties. At 
frequencies > 100 GHz, thermal dust emission dominates and is 
well described by a modified blackbody with an emissivity index 
/5 = 1.71 ± 0.12 and dust temperature = 25.7 + 1.7 K. 

We note that the lack of AME in M42 is likely to be due 
to the diff'erence in conditions (e.g., intense radiation field) that 
may result in depletion of the smallest dust grains, which would 
result in a large reduction in the amount of spinning dust emis- 
sion. This is a diff'erent situation to the more diffuse warm 
ionised medium (WIM) studied by iDobler et"aD (l2009h , where 
they detected AME in the WIM at the level of ^ 20 % relative to 
the free-free emission. 

Eig. [8] shows the spectrum of AME-G353.03+16.90 in the 
p Oph West region; the flux densities are listed in Table O The 
same model used for the Perseus region (Sect. l4.2l) was fitted to 
the data, except for the details of the spinning dust component, 
which will be discussed in Sect.[6j Significant CO line contami- 
nation is visible at 100 and 217 GHz, so these two bands were not 
included in the fit. There is minimal free-free emission within the 
aperture, and hence the error at low (0.4-2.3 GHz) frequencies 
is dominated by fluctuations in the local background, as indi- 
cated by the large error in the 2.3 GHz flux density. Note that the 
408/1420 MHz data were included in the fit but were consistent 
with zero. The best-fitting parameters are: EM =1+48 cm"^ pc; 
Td = 20.7 + 1.6 K; = 1.75 + 0.18; and CMB contribution 
ArcMB = 82 + 61 |LtK. The fit of the overall model is very good 
(x^/dof = 0.29). 

As a test of the impact that the Planck data was having on 
the spectrum, we repeated the analysis omitting the Planck data. 
As expected, the impact was considerable. In particular, the ther- 
mal dust model is constrained much more strongly the HEI data 
are included. Eor example, the uncertainty in the dust tempera- 
ture increases to +2.5 K, while the uncertainty in the emissivity 
index increases to +0.36. This is crucial for making precise mea- 
surements of the AME at frequencies, particularly above 50 GHz 
where the subtraction of thermal dust is critical. 

In Eig. [9] we show the residual spectrum after subtracting the 
free-free, CMB, and thermal dust models. The additional error 
from the subtraction of the models is added in quadrature to the 
flux density errors to obtain the uncertainty on the residual flux 
densities, listed in Table [3l As in Perseus, the residual spectrum 
has a clearly defined convex shape, peaking at a frequency of 
30 GHz; the excess is significant at the 8.4(T level. The spec- 
trum at 50-100 GHz is flatter than a single spinning dust model 
can easily account for, suggesting a different environment or dis- 
tribution of dust grains. This may also be due to variations in 
the dust emissivity index due to multiple dust components (e.g., 
Einkbeiner et al. 1999). To test this possibility would require the 
inclusion of the 100/217 GHz data after careful subtraction of 
the CO line contamination. 
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Fig. 6. Maps of the p Ophiuchi molecular cloud region at their original angular resolution. From left to right, top row: Planck 28.5; 
44.1; 70.3; and lOOGHz. Bottom row: Planck 143 and 857 GHz; 1.4GHz; and Ha. The maps cover 5° x 5° centred on (/,/?) = 
(353?05, +16?90) and have linear colour scales. The graticule has 1° spacing in Galactic coordinates. The circular aperture and 
background annulus, which were used to calculate the flux density (see text), are indicated. The strong AME is evident at 28.5 and 
44.1GHz. 
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Fig. 8. Spectrum of AME-G353. 05 +16.90 in the p Ophiuchi 
West molecular cloud. The best-fitting model consisting of 
free-free (orange dashed line), spinning dust, CMB (black 
double/triple-dotted line), and thermal dust (light blue dashed 
line), is shown. The spinning dust model consists of two compo- 
nents: high density molecular gas (magenta dot-dashed line); and 
low density atomic gas (green dotted line). The 100/217 GHz 
data are contaminated by CO line emission and are not included 
in the fit. 
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Fig. 9. Spectrum of AME-G353. 05 -h 16.90 in the p Ophiuchi 
molecular cloud after subtracting the best-fit free-free, CMB 
and thermal dust components. A theoretical spinning dust model 
consisting of two components is shown as the black solid line; 
the magenta dot-dashed line is for dense molecular gas and the 
green dotted line is for irradiated low density atomic gas (see 
Sect.©. The 100/217 GHz data are contaminated by CO line 
emission and are not included in the fit. 



A theoretical spinning dust model consisting of two com- 
ponents (associated with atomic and molecular gas) is over- 
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Table 3. Flux densities for AME-G353.05+ 16.90 in the p 
Ophiuchi molecular cloud and residual flux densities when free- 
free, CMB and thermal dust components are removed. 



Frequency 


Flux density 


Flux density residual 


[GHz] 


[Jy] 


[Jy] 


0.4 


-6.8 + 9.4 


-7 + 11 


1.4 


-1.0 + 6.8 


-1.0 + 8.3 


2.3 


1.3 + 6.9 


1.2 + 8.3 


22.7 


26.3 + 5.5 


24.8 + 6.6 


28.5 


29.6 + 5.6 


27.3 + 6.6 


33.0 


30.7 + 5.3 


27.2 + 6.3 


40.7 


27.7 + 4.6 


21.9 + 5.8 


44.1 


27.0 + 4.4 


19.9 + 5.7 


60.6 


26.3 + 4.5 


9.8 + 6.5 


70.3 


31.4 + 5.1 


6.5 + 7.8 


93.4 


63.6 + 8.9 


6+15 


100 


84 + 13 


13 + 19 


143 


222 + 31 


5 + 49 


217 


989 + 16 


130 + 200 


353 


4100 + 670 


-180 + 780 


545 


17100 + 2900 


400 + 3200 


857 


58300 + 8800 


500 + 10000 


1249 


122000 + 22000 


-5000 + 25000 


2141 


218000 + 37000 


20000 + 41000 


2997 


125000 + 26000 


-9000 + 32000 



plotted and provides an excellent fit to the data. Note that the 
denser molecular gas represents the dominant AIME with peak at 
^ 30 GHz, while the irradiated low density atomic gas accounts 
for the 50-100 GHz part of the spectrum. The physical model is 
presented in Sect.[6l 

6. Modelling the anomalous microwave emission 
with spinning dust 

Compar isons of the AlVIE with dust emission in localised re- 
gions (Casassus et al. "20061: IScaife et al.ll2009l: lTibbs et al." 20101: 
Castellanos et a l. 2011) and in all-sky surveys (Lagache 20031; 
Ysard et al.l I2Q1Q ) suggest small dust grains (polycyclic aro- 
matic hydrocarbons, PAH) as carrie rs of this emis sion, sup- 
porting the spinning dust explanation (iDraine & Laza rian 1998). 
The present and previous measurements of the AME show 
that this component peaks at 20-40 GHz with a brightness of 
0.01 to 0.1 MJy s r~^ (see references in Sect.[I])- Current spin- 
ning dust model s ("Ali-Hai'moud et al."2009'; lYsard & Verstraetd 
I201QI: iHoang et al. 2010; Silsbee et al. 201 1^) indicate that low- 
density (^H < 30 cm"^), atomic diff'use gas peaks in this fre- 
quency range. Therefore the AME has often been explained 
with spinning dust emission arising mostly from regions of low 
density gas, such as from the Cold Neutral Medium (CNM) or 
Warm Neutral Medium (WNM). The column density A^h asso- 
ciated with the AME can be derived from observations provided 
the emissivity of spinning dust is known: adopting an average 
of 10"^-^ MJy sr"^ cm^ H"^ from theoretical models, we find 
that A^H ~ 10^^ cm"^ for a flux density level of 0.1 MJy sr"^ For 
a gas density of 30 cm"^, the extent of the emitting region is 
thus Nu/nu -100 pc, much larger than the observed size of the 
present AME clouds (< 10 pc). The situation is even worse for 
the warm ionised medium (WIM) component (^h < 1 cm"^). 
The spinning dust emission of these regions must therefore in- 
clude a significant contribution from g as which is denser by 
at l east a factor of 10, as s uggested by ICasassus et al.l (120081) 
and ICastellanos et aP (l201lh . In the preliminary models of the 



Table 4. Parameters in the spinning dust models for Perseus and 
p Ophiuchi: xu and xq are the abundances of H^ and C^ ; j is the 
molecular fraction 2^h2/^h; (^o is the centroid of the PAH size 
distribution; be is the PAH abundance of carbon in PAHs; and z 
is the depth of the emiss ion region al o ng the line of sight. The 
radiation field is that of iMathis et al.l (1 19831) multiplied by the 
factor Go. The ionised gas accounts for the free-free emission 
but does not contribute to the spinning dust emission. 



Gas state 


Molecular 


Atomic 


Ionised 


Perseus 


Nu cm-'] 


11.7 


1.3 


0.4 




250 


30 


1 


z [pc] 


15.1 


14.0 




Go 


1 


2 




r[K] 


40 


100 


8x 10^ 


xu [ppm] 


112 


410 


10^ 


xc [ppm] 


<1 


100 




y 


1 


0.1 




ao [nm] 


0.58 


0.53 




be [ppm] 


68 


68 








1.65 








18.5 




T250 




9.4 X 10-4 




p Ophiuchi 


A^H [10^^ cm"^] 


18.2 


0.4 


0.4 


nn [cm-^] 


2x 10^ 


200 


0.5 


z [pc] 


0.3 


0.6 




Go 


0.4 


400 




T[K] 


20 


103 


8x 10^ 


xu [ppm] 


9.2 


373 


10^ 


Xc [ppm] 


<1 


100 




y 


1 


0.1 




ao [nm] 


0.60 


0.38 




be [ppm] 


65 


50 








1.75 




rd[K] 




20.7 




^250 




3.2 X 10-3 





Perseus and p Ophiuchi regions presented below, we show that, 
assuming plausible physical conditions and PAH size distribu- 
tions, spinning dust emission from dense gas can explain most 
of the AME. 

The rotation of PAHs in spac e is governed by IR emis - 
sion and gas-grain interactions dPraine & LazarianI 1 19981) . 
Quantitative modelling of spinning dust emission requires the 
knowledge of many (~10) parameters describing the state of the 
gas, the radiation field, and the PAH size distribution , as well 
a s the electr ic dipole moment. A s shown by Ali-Haim oud et al.l 
(I2QQ9I) and lYsard et al.1 (l201Ql) . the peak frequency of spin- 
ning dust emission in neutral gas (and for radiation intensity 
Go < 10^, where Go is the multiplying factor relative to the 
mean rad iation intensity in the interstellar medium, as defined 
by M athis et al.llT983) depends chiefly on the gas density and on 
the size of PAHs0 The largest PAHs have the lowest emissivity. 
In low temperature gas (Go < 3 and > 30 cm"^), the gas- 
grain collisions dominate the excitation and the damping of the 
grain rotation. Thus the spinning dust spectrum becomes sen- 
sitive to the abundance of the major ions, H^ and C^, denoted 
xu and xq respectively. To compute the spinning dust emissivity. 



^ We assume here t he electric dipole moment to be as in 
iDraine & LazarianI 1 19981 . a prescription al so shown to be c ompatible 
with the AME extracted from WMAP data (lYsard et al.ll201Qr) . 
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we use SPDUSlH in the case 2 of llSilsbee et alJ (l201lh . Rather 
than fitting all the spinning dust parameters on the observed 
spectrum, we derive realistic values from observations, and fit 
the amplitude to the spectrum. We derive xh from the ionisation 
balance of hydrogen assumi ng a standard cosmic-ray ionisation 
rate of 5 x 10"^^ s"^ H'^ (IWilliams et al.lll998l: IWolfire et al.1 
1200 3). Conversely, xq cannot be estimated so readily because 
C undergoes reacti ons with H2 and PAHs (Rolligetal. 20061; 
IWolfire etal.ll2008h : we therefore take xq as a free parameter. 
The size distribution of PAHs is assumed to be a log-normal of 
centroid and width a = 0.4. The abundance of PAHs is set by 
be, the abundance of carbon locked up in PAHs. 

Now we describe the physical conditions within our two re- 
gions, the Perseus and p Ophiuchi molecular clouds, as derived 
from observations, and how we take them into account in our 
modelling (see Table|4l). For both regions, the free-free contribu- 
tion is that of the H 11 region and the thermal dust emission is de- 
scribed as a modified blackbody, ly = T250 (/l/250|Ltm)"^ By(T) 
where T250 is the dust opacity at 250 jam, T the dust temper- 
ature (see Table |4l), and By is the blackbody brightness. We 
do not include the contribution of the ionised gas to the spin- 
ning dust emission because of its low density and also because 
PAHs have been shown to be destroyed efficiently in such envi- 
ronments (Lebouteiller et al. 2007, 2011). In addition, we note 
that the current SPDUST code does not include ionising photons 
(hv > 13.6 eV), which are important for the PAH charge and IR 
emission. 



6.1. Perseus molecular cloud 

At the degree scale, the Perseus line of sight intersects a 
molecular cloud surround ing the Hii region associated with 
the HD 278942 BOV star. iRidge et al.l (|2006|) argued that this 
ionised shell has a density of ~ 1 cm"^. The gas density and 
temperature of the molecular component are taken from the ana l- 
ysis of observed C2 lines carried out byE slesias- GrothI (I2QTTP I 
The parameters of our model, which provides a good match to 
the AME, are given in Table ID and the corresponding emission 
is shown in Fig. \5\ We note that the depth, z, of the emissive 
region along t he sightline is coin patible with determinations at 
smaller scales (iRidge et al.ll2006l) . The spinning dust emission is 
dominated by the contribution of the dense gas; the low-density 
atomic gas appears to play only a minor role. Thi s is consistent 
with the results of iPlanck Collaboration! (1201 Iqb . which show 
that, on large scales, there is bright AME associated with molec- 
ular clouds traced by CO everywhere in the Galaxy, while the 
association with the diffuse and atomic phase correlated with H i 
is less clear. 

In contrast to previous work where xu = was assumed in 
molecular gas, our ionisation balance of H gives xh ~ 3 x 10"^, 
a change particularly important for Perseus. We take xq to 
be 1 ppm: at this level the C^ ion has no inffuence on the 
spinning dust emission. Adopting larger sizes for the PAHs 
(ao = 0.58 nm) than those used in previous spinning dust mod- 
els [but comparable to the diffuse interstellar medium (0.64 nm); 
ICompiegne et al.|[20TL] . the molecular gas accounts well for the 
AME. 



6.2. The p Ophiuchi molecular cloud 

The p Ophiuchi spectrum includes contributions from a photo- 
dissociated interface region (PDR) and the ionised gas associ- 
ated with the star HD 147889. We note that IR observations have 
shown that the PAH emission arises from dens e irradiated inter- 
fa ces (lAbergel et al ."2002l: iHabart e"taDl2005D . As constrained 
by lHabart et al.l (2003), the PDR has a density gradient between 
the star and the molecular cloud. For the present preliminary 
model, we represent this interface with two effective media at 
low (200 cm"^) and high (2 x 10^ cm"^) densities irradiated by 
strong and weak radiation fields. In the more excited p Ophiuchi 
region we are able to match the AME by requiring: (a) an at- 
tenuated radiation field in the dense gas; and (b) smaller and 
less abundant PAHs in stronger radiation fields, in line with re- 
cent work suggesting t hat PAHs are pho t oevaporated froni larger 
carbonaceous grain s (Bern e et ani2007l: ICompiegne et al ] l20Q8l: 
IVelusamy & Lang er 2008). The PAH abun dances we assume fo r 
both regions are within the range found by H abart et al.l (l2003h . 
The fit is shown in Fig. [9] and the parameters of our model are 
summarised in Table IH 

We note, however, that the determination of ao. Go and nu 
is degenerate in the present modelling; you can create almost 
identical spinning dust curves by varyi ng these p arameters in 
different ways (see, e.g., Al i-Haimoud et al .1120091) . Lifting this 
degeneracy involves treating the radiative transfer along with the 
spinning dust motion. The present results are therefore sugges- 
tive only and we emphasise that the physical conditions of our 
models are derived from scales much smaller than the resolu- 
tion of our observations. The previous discussion also indicates 
that quantitative model spectra can only be obtained from a con- 
sistent treatment of the gas state, the radiative transfer, and the 
spinning motion of dust grains. Finally, a combination of multi- 
frequency data, including detailed IR measurements, is clearly 
needed to extract detailed information about the environment 
and the dust grains. 



7. New regions of anomalous emission 

We have searched the Planck maps for new regions of AME on 
the scale of ~ 3° or smaller. To do this, one must suppress the 
other emission components that contribute to the frequencies 
where AME is the strongest (~20-60GHz). Numerous meth- 
ods can be used (e.g.. Leach et al. 2008; Delabrouille & Cardoso| 
2009). However, we have taken a simplistic approach to identify 
bright AME regions. We subtracted a simple spatial model of 
known emission mechanisms (synchrotron, free-free, and ther- 
mal dust) using extrapolations of existing templates from obser- 
vations or theoretical expectations, and then examined the loca- 
tions of the remaining emission in more det ail. We used map s 
smoothed to 1°, with HEALPix A^side = 256 (Gors ki et al]l2005h . 

We subtracted the synchrotron emission fro m the CMB- 
subtracted Planck maps using the 408 M Hz map ( Haslam et al.J 
1982 ), while the 1 420 MHz ( Reich 1 982t iReich & Reich 1986: 
Reic h etal.ll2QQl1) and the 2326MHr (|jonas et al.lll998h maps 
were used in conjunction with 408 MHz data to determine the 
spectral index. Brightness temperature spectral indices were re- 
stricted to lie in the range -2.0 > p > -4.0 (T oc y^). For 
the gap region where no 1420/2326 MHz data exist, we assume 



P = -3.0. We used ve rsions of the maps that have had no 

4 http : //www . tapir . caltech . edu/~ yacine/spdust/spdust . htifl^>int sources subtracted (IPlatania et al]|2003|) so as to allow di- 

5 A single high resolution measurement may not be representative rect comparison with t he Planck data. For fre e-free emission, 
of the line-of-sight on degree scales, since the density and temperature we use the H a map of iDickinson et al.1 (l2003l) . which is based 
may vary considerably across the beam. on WHAM (iHaffner et al.ll2003l) and SHASSA (iGaustad et al.l 
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Fig. 10. Residuals in the full sky Planck LFI 28.5 GHz 1° smoothed map after subtraction of synchrotron, free-free and thermal dust 
emission (see text). 12?5 x 12?5 cut out maps are shown for the Perseus and p Ophiuchi molecular clouds, and the two new regions 
of AME, AME-G107. 1+5.2 and AME-G173.6+2.8. The graticule spacing is 5° in Galactic coordinates. 



|2QQ1|) data corrected for dust absorption. The map is also used 
to correct the synchrotron maps for free-free emission, assuming 
an average electron temperature of 8000 K. The Rayleigh- Jeans 
tail of the thermal (vibrational) dust component is determined 
from the Planck maps at 143 and 545 GHz, thereby avoiding 
significant CO contamination. The 545 GHz data were used for 
normalisation and the 143 GHz data were used for the spectral 
index of the power-law. This was constrained to lie between 
and +3 to avoid extrapolation artefacts in some areas at high lat- 
itudes where the uncertainty on the spectral index is large. 

The full- sky map of the residuals from the 28.5 GHz 
Planck LEI data is shown in Eig. [TOl Artefacts from the sub- 
traction are present. The division between the 1420 and the 
2326 MHz maps is evident as well as the lack of observations 
around the south celestial pole. Other large-scale features are 
also present. These are extended features and hence do not aff'ect 
our search for compact areas of AME. A large number of AME 
candidates are evident in the map, including the strong Perseus 
and p Ophiuchi AME regions described in Sects. IH and [3 We 
examined ~ 50 regions in detail. 

The next step in identifying new AME regions is to con- 
struct a spectrum for each candidate, using aperture photometry. 
In each case there is strong thermal dust and free-free emission 
indicative of an associated H ii region. We now illustrate the re- 
sults of our search in two of the compact regions which show 
excess emission in the ~20-60 GHz range. 



7.1. AME-G1 73.6+2.8 

AME-G173.6+2.8 (05Ml^, +35°5r J2000.0) is a strong, com- 
pact region of dust e mission to\yards t he Galactic anticentre. 
The Hii region S235 (ISharplessI 11959^ lies within this area. 
It is a region of massive star formation with a young stellar 
object at its centre, and it has several compact radio compo - 
nents within several arcmin (iNordh et al.ll984l:lEelli et al.l2006h . 
Radio recombination li ne observations indic ate a radial velocity 
of -24.5 ± 0.7 km s"^ (Paladini et al]l2003l) . S235 lies at a dis- 
tance of 1.8 kpc and is therefore within the Perseus spiral arm. 

Maps of the region at a range of frequencies are shown at 
the observing resolution of each in Eig. [TTJ The AME is strong 
at 28.5 GHz while the dust distribution at the higher Planck fre- 
quencies shows the small-scale structure of the dust within the 
28 .5 GHz beam. It c an be seen that the Ha at the 6' resolution 
of lEinkbeined (l2003l) is off'set from the dust, most likely due to 
absorption by the dust. 

The spectrum of the 1?3 degree (COSMOSOMAS) region 
centred on AME-G173.6+ 2.8 is shown in Eig.fHl The radio flux 
densities are at 820 MHz ( Berkhuiisenlll972h . 1.4 GHz (ReichI 
Il982) and 11-17 GHz (COSMOSOMAS). The spectrum and 
residuals of the region are derived using the COSMOSOMAS 
filtering method described in Sect. l4.2[ 



The spectrum is well-fitted using a combination of thermal 
dust (with r = 19.6±0.6Kandanemissivity indexofyS = 1.54 + 
0.07), free-free emission (EM = (2.4 7 + 13) x lO^cm'^pc) 
and spinning dust (Dra ine & LazarianI (1 19981) CNM model with 
A^H = (2.7 + 0.3) X 10^^ cm"^), with a smaU CMB contribution of 
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Fig. 11. Maps of the AME-G173.6+2.8 region: 28.5 GHz (linear 
scale); 70.3 GHz (linear scale); 545 GHz (logarithmic scale); and 
Ha (logarithmic scale). The maps are 8° x 8° and the graticule 
spacing is 1° in Galactic coordinates. 



ArcMB = 130 + 40 |LtK. The residual spectrum shows the signa- 
ture of spinning dust peaking in the frequency range 20-30 GHz, 
detected with a significance level of 6 Act. The main uncertainty 
in the AME spectrum is the level of free-free emission to be sub- 
tracted. As shown, the spinning dust emission in the 20-30 GHz 
region is ~ 2 times the free-free emission; the uncertainty is 
< 10-20 %. The current spectral coverage defines the peak of the 
spectrum and provides useful constraints in modelling the AME. 
An ultracompact H ii region, producing free-free emission that is 
optically thick at frequencies < 10 GHz, may be contributing to 
the flux density at > 10 GHz. However, it cannot easily account 
for the shape at high (> 40 GHz) frequencies. 

7.2. AME-G1 07.1+5.2 

AME-G107.1+5.2 (22^22^, +63°23^ J2000.0) is a dust complex 
containing the star-forming Hii region S140 (I SharplesslfT959h 
with Vlsr = -8. 5 + 1.0 km s~^ Th e region lies at the edge of the 
Cepheus bubble (iPatel et al.l ll998) at a distance of 800 pc. 

Maps of the region are shown in Fig.[T3]using the same for- 
mat as in Fig. [TT] Dust is widespread over the region but with 
a concentration at the position of the AME. There are multiple 
dust components within the 28.5 GHz beam. There is a small ofl'- 
set of the 28.5 GHz feature of 15' to lower latitudes and higher 
longitudes than the brightest feature in the dust maps. This indi- 
cates a higher radio emissivity of the dust in this ofl'set position. 
There is also an ofl'set in the Ha emission, again possibly due to 
absorption by the associated dust. 

The spectrum estimated for AME-G107. 1 +5.2 and the resid- 
ual spectrum are shown in Fig. [141 The least- squares fit to the 
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Fig. 12. The spectrum (top) and residuals (bottom) for AME- 
G173.6-h2.8 after removing free-free emission (orange dashed 
line), thermal dust emission (light blue dashed line) and CMB 
anisotropics (not visible). A spinning dust model is shown as a 
magenta dot-dashed line. 



free-free emission (EM = 150 + 39 cm ^ pc), thermal dust (Td = 
19.0± 2.6 K, /5 = 2.04 ± 0.45), spinning dust toraine & LazarianI 
(1 19981) WIM model with A^h = (1.76 ± 0.44) x lO^^ cm'^), 
and CMB (AJcmb = 75 ± 38 ^iK) is well-defined. The AME 
is detected at a significance of 4.8(T with a shape that is con- 
sistent with spinning dust. However, an ultracompact H ii region 
producing optically thick free-free at frequencies below 10 GHz 
may also be contributing. 

7.3. Assessment of the search of new AME regions 

We have shown that it is possible to detect bright compact AME 
regions in the Planck data using a rather simplistic selection pro- 
cess. In order to verify the identifications it is necessary to add 
ancillary data at radio and FIR frequencies. Now that identifica- 
tions have been made, new ancillary data covering a larger fre- 
quency range will be sought. Further development of the search 
algorithms will also be worthwhile. 

The strength of AME relative to the free-free and thermal 
dust emission in these two regio ns can be co mpared with the 
values for Hii regions obtained by lTodorovic et al. (2010) in the 
VSA survey of the / = 27°-47° section of the Galactic plane. 
The latter survey found the ratio of 33 GHz AME to free-free 
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Fig. 13. Maps of the AME-G107. 1+5.2 region: 28.5 GHz (linear 
scale); 70.3 GHz (linear scale); 545 GHz (logarithmic scale); and 
Ha (logarithmic scale). The maps are 8° x 8° and the graticule 
spacing is 1° in Galactic coordinates. 



flux density to be 0.40 ± 0.11; the corresponding values for 
AME-G173.6+2.8 and AME-G107. 1+5.2 are 0.67 and 1.4, re- 
spectively. In the case of the ratio of 33 GHz AME to 100 |Ltm 
thermal dust emission, the VSA survey value was 1.1 x 10"^ 
(or 2 X 10""^ at the 15 GHz peak of the spectrum), while the 
AME-G173.6+2.8 and the AME-G107. 1+5.2 values are some- 
what larger at 3.3 x 10"^ and 3.8 x 10"^, respectively. It is 
however not surprising that the 30 GHz AME is stronger in the 
two Planck H ii regions than in the VSA H ii regions, because 
the former were selected as the brightest AME sources in the 
Planck maps. The contribu tion from ultracompact H ii regions 
(IWood & Churchwelll[l989h may also be a significant contribu- 
tor to the AME at low Galactic latitudes. 

Finally, we note that the derived dust temperature for these 
Hii regions, at Td ^ 19 K, is closer to that of the diff'use inter- 
stellar medium than to the warm d ust (Td ~ 30-50 K) usuall y 
found surrounding OB stars (e.g., IWood & Chur chwell 1989). 
The dust observed along the line of sight at ~ 1° resolution is 
likely to be a mixture of populations, including higher density 
regions (e.g., molecular clouds) that could dominate the signal 
on these scales. This may explain the lower average temperature 
the presence of significant AME in these sightlines. 

As selection procedures develop, we expect many AME re- 
gions to be discovered in a range of physical environments, 
which will provide data necessary to acquire a full understanding 
of the AME emission mechanism. 
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Fig. 14. The spectrum (top) and residuals (bottom) for AME- 
G107.1-h5.2 after removing free-free (orange dashed line), ther- 
mal dust emission (light blue dashed line) and CMB anisotropics 
(black double-dotted line). A spinning dust model is shown as a 
magenta dot-dashed line. 



8. Conclusions 

Filling in the gap between earlier radio and FIR measurements, 
the frequency coverage of Planck has provided a unique oppor- 
tunity to establish comprehensive spectra of regions of AME. 
Planck has revealed the high frequency side of the spectral peak. 
The new spectra are the basis for understanding the emission 
mechanism and the environment in which it occurs. The evi- 
dence from the present observations strongly favours the spin- 
ning dust mechanism (electric dipole radiation). Planck provides 
a rich data set that can be used as a basis for developing a realis- 
tic understanding of the AME mechanism in a range of Galactic 
environments. 

The two best- studied AME sources that have extensive an- 
cillary data are the Perseus and p Ophiuchi molecular clouds. 
In each, the spectrum is well fitted by free-free, thermal dust, 
and spinning dust, with a small contribution from the CMB. 
Spinning dust provides a good fit to the microwave (10- 
100 GHz) part of the spectrum, which peaks at ^ 30 GHz. 

Theoretical spinning dust spectra are presented for a physical 
model consisting of molecular and atomic states. It is possible to 
derive physical parameters that are consistent with the environ- 
ment and still provide a good fit to the data. Using parameters 



14 



Planck Collaboration: Planck Early Results. XX. New light on anomalous emission from spinning dust grains 



constrained at smaller angular scales, the 20-40 GHz AME peak 
in Perseus is well explained with spinning dust emission aris- 
ing from dense, molecular gas (^h > 200 cm"^) subjected to a 
few times the typical interstellar radiation field. Low-density gas 
is only a minor contributor to AME in Perseus; however, in p 
Ophiuchi, although dense gas accounts for the peak at ^ 30 GHz, 
irradiated low-density atomic gas may be contributing in the 
range 50-100 GHz. The picture seems to be that smaller PAHs 
are found in PDRs (Go > 100), as suggested by recent Spitzer 
observations. However, the determination of the PAH size is de- 
generate with that of and Go and quantitative conclusions will 
only be obtained from consistent modelling of the gas state, ra- 
diative transfer, and spinning dust. At this level of modelling it 
is not possible to constrain the electric dipole moment of PAHs. 
Future work to carry out more detailed modelling is required. 

Our preliminary search for new AME regions in the 
Planck data has revealed many candidates. They were uncov- 
ered by subtracting synchrotron, free-free and thermal dust emis- 
sion based on the usual spatial templates. Two of the new candi- 
date regions that show AME at > 5 cr significance have spectra 
consistent with a spinning dust spectral shape. Additional multi- 
frequency data (e.g., s uch as the 5 GHz C-Band All- Sky Survey, 
C-BASSQ; i King etal.l[2QT 0) as well as high-resolution observa- 
tions will be needed to understand the detailed structure of the 
AME in these regions. 
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